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ABSTRACT 

We highlight the impact of cluster-mass-dependent evolutionary rates upon the evolution of 
the cluster mass function during violent relaxation, that is, while clusters dynamically respond 
to the expulsion of their residual star-forming gas. Mass-dependent evolutionary rates arise 
when the mean volume density of cluster-forming regions is mass-dependent. In that case, 
even if the initial conditions are such that the cluster mass function at the end of violent 
relaxation has the same shape as the embedded-cluster mass function (i.e. infant weight-loss 
is mass-independent), the shape of the cluster mass function does change transiently during 
violent relaxation. In contrast, for cluster-forming regions of constant mean volume density, 
the cluster mass function shape is preserved all through violent relaxation since all clusters 
then evolve at the same mass-independent rate. 

On the scale of individual clusters, we model the evolution of the ratio between the dy- 
namical mass and luminous mass of a cluster after gas expulsion. Specifically, we map the 
radial dependence of the time-scale for a star cluster to return to equilibrium. We stress that 
fields-of-view a few pc in size only, typical of compact clusters with rapid evolutionary rates, 
are likely to reveal cluster regions which have returned to equilibrium even if the cluster expe- 
rienced a major gas expulsion episode a few Myr earlier. We provide models with the aperture 
and time expressed in units of the initial half-mass radius and initial crossing-time, respec- 
tively, so that our results can be applied to clusters with initial densities, sizes, and apertures 
different from ours. 
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1 INTRODUCTION 

Star cluster formation conditions have a profound impact on the 
properties of young star cluster systems. In particular, after mas- 
sive star activity has cleansed a cluster-forming region (CFRg) of 
its residual star-forming gas, the mass fraction of stars bound to 
a cluster decreases steadily with time (infant weight-loss) until the 
cluster either reaches a new equilibrium, or is completely disrupted. 
This dynamical response of a cluster to gas ex pulsion is called vi- 
olent relaxation. iBaumgardt & Kroupa I d2007l ) express the bound 
mass fraction of stars, F^, as a function of the CFRg star formation 
efficiency at the onset of gas expulsion, SFE, of the gas expulsion 
time-scale, T geX p/t cross , and of the impact of an external tidal field 
which they quantify through the ratio of the initial half-mass radius 
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to the initial tidal radius of the embedded-cluster, r/Jr,. That is: 

Ffr — Ffoit J ' T cross ,SF E , Tgexpl ^crossi ^h/^t^) ? (1) 

where Fj, is denned as the mass fraction of stars within the 
cluster tidal radius at time t. Note that the time-scale for gas 
expulsion, Tg CTp , and the time since gas expulsion, t, are both 
expressed in units of the CFRg crossing-time, Z C ross- This im- 
mediately implies that denser CFRgs give rise to faster-evolving 
-1/2 

clusters since T crosj 0= P CFRg , with PcFRg the CFRg mean volume 
density. This property is the crux of the present contribution. 
The higher the SFE, the longer the gas-expulsion time-scale 
tgexpl Across, the weaker the tidal field impact r/Jr, (i.e. the 
deeper the cluster lies within its limiting tidal radius), the smaller 
the cluster infant-weight loss, 1 — F/,, at time t/t cross . Cluster 
violent relaxation therefore leaves an imprint on the age distri- 



bution l lParmentier & Fritze" 20091) and mass function of young 
star clusters dKroupa & Boilv I 120021: IBaumgardt et alj [20081 : 



IParmentier et al . 2008; Pa rmentier & Kroupa 11201 JT 
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During the few 10s of Myr after gas expulsion, infant weight- 
loss drives the cluster mass function towards lower cluster masses. 
The amplitude of this leftward shift as a function of time, and mass 
(i.e. does the cluster mass function shape vary with time?), consti- 
tutes a powerful probe into cluster formation conditions. Most ob- 
servational evidence gathered so far seem to support the scenario 
of mass-independent infant weight-loss, 1 — Fj, at that young age. 
When plotted as the number of clusters per linear mass interval, the 
young cluster mass function in the present-day Universe is a power- 
law of s lope ~ —2, dN r i qc m^ 2 dm r /, irrespective of the cluster 
age (e.g.lLada & Lada Il2003l ; lOev et alj|200l iDowell et alj|2008l ; 
IChandar et al.l201(I) . It should be kept in mind that, in assessing the 
slope of power-law cluster mass functions, systematic errors (e.g. 
the errors inherent to the derivations of cluster mass estimates by 
different research groups) are significantly larger t han the statistical 
ones (i.e . the formal fitting errors). For instance, iBaumgardt et al.l 
din prep. I) find that the mass function slope of the Large Magel- 
lanic Cloud clusters younger than 200 Myr is —2.32 ±0.11. T his 
is steeper than the slope reported bv lde Griis & Anders I d2006l) for 
a similar cluster age range, i.e. —1.94 ±0.10. This simple exam- 
ple shows that to detect an intrinsic change of ~ 0.2 in the cluster 
mass function slope is not straightforward. We note that through 
a careful analysis of the combined star cluster systems of several 
galaxies - an approach whic h aims at decre asing statistical uncer- 
tainties at the high mass end -[Earsen ( 2009) shows that the cluster 
mass function obeys a Schechter function, rather than a featureless 
power-law, with a cut-off mass higher in starbursts and mergers 
than in quiescent spiral discs. 

The amount of mass-independent dissolution affecting young 
clusters remains heavily debated. Based on the relation between 
the mass of the most massive cluster and the age range sam- 
pled, iGieles & Bastian ] d2008h infer almost negligeable mass- 
independent di ssolution, i.e. 20 % pe r age dex as compared to 90 % 
per age dex for lChandar et all d2010h . We stress that these observa- 
tional uncertainties do not affect the general conclusions we present 
here. 

If infant weight-loss is mass-independent, this straightfor- 
wardly implies mass-independent SFE and mass-independent 
tgexp/tcmss, although the const raint on T eexn I Across is loose r than 
for the SFE (see section 4.1 in lParmentier & Kroupa ||201 ll for a 
discussion). For a given external tidal field (e.g. a limited region of a 
given galaxy), it also implies the mass-independence of r/,/r t . This 
constraint is robustly satisfied if the mass-ra dius relation of CFRgs 
is one of constant mean volume density JParmentier & Kroupa I 
l201ll) . i.e. rcFRg 01 m CFRg' w ^ m r CFRg and mcFRg the radius and 
mass of CFRgs, respectively^ Mass-independent SFE, % ge xp/ Across 
and Th/rt therefore imply that the cluster mass functions at the on- 
set and end of violent relaxation have the same shape, albeit differ- 
ent amplitudes due to cluster infant weight-loss. 

A fourth, and so far overlooked, aspect is how the cluster 
evolutionary rate shapes the cluster mass function during violent 
relaxation. If the mass-radius relation of CFRgs is one of con- 

1/2 

stant mean surface density (i.e. r^pR? <* m CFRg>' more massive 
CFRgs have a lower mean volume density and, therefore, evolve 
on a slower time-scale. Conversely, if CFRgs all have the same ra- 
dius, r^pRg, irrespective of mass, higher-mass CFRgs give rise to 



faster- evolving star clusters, as illustrated by fig. 1 of IParmentierl 
In that example, the constant radius also implies that more 
massive CFRgs have deeper gravitational potential wells hence 
slower gas-expulsion time-scales. This allows high-mass clusters 
to survive despite a low SFE of 0.20, while low-mass ones are 
disrupted. The mass-dependent gas-expulsion time-scale thereby 
carves a turnover in the initially power-law cluster mass function. 
The same figure also shows that, as clusters evolve after gas ex- 
pulsion, transient patterns sets in the cluster mass function due to 
the mass-dependent evolutionary rate. At an intermediate age of 
15 Myr, clusters less massive than 3000 Mq (i.e. low-density clus- 
ters) have hardly evolved and their mass function sticks to the initial 
one. That is, none of their stars have yet crossed their tidal bound- 
ary. Conversely, clusters more massive than 10 4 Mq (i.e. high den- 
sity clusters) have completed their violent relaxation, and the bell- 
shaped cluster mass function already shows up in the high-mass 

regime. 

In this contribution, we expand the work of IParmentierl 

( 120101) . We highlight how the CFRg mass-radius relation affects 
the evolutionary rate of star clusters in dependence of their mass 
and, therefore, the cluster mass function evolution. In particular, we 
quantify transient mass-dependent effects which do take place even 
if the cluster mass function shapes at the onset and at the end of 
violent relaxation are alike. Although we illustrate our point with a 
power-law cluster mass function experiencing 90 % infant weight- 
loss, we emphasize that the cluster mass function patterns high- 
lighted below do take place regardless of the initial cluster mass 
function shape and regardless of the infant weight-loss rate. 

While the cluster mass function evolution constitutes a good 
tracer of violent relaxation on the scale of star cluster systems, the 
dynamical evolution of an individual cluster after gas expulsion can 
be traced by the ratio of its dynamical-to-luminous mass estimates, 
Mdxn/Mi um . Again, the shorter the initial crossing-time, the faster 
the return of the cluster to virial equilibrium where M^ yn jM[ um ~ 1 . 
The question of whethe r a young cluster can be observed in virial 
equilibrium (e.g. Wd 1, iMengel & Tacconi-Garman |2007|) despite 
having been put out of equilibrium by a significant gas-expulsion 
episode in the recent past is thus directly related to its initial 
crossing-time hence to the CFRg mass-radius relation. As we shall 
see in the present work, the problem is complexified by the size of 
the aperture with which the cluster is observed. For a given time- 
span after gas expulsion, smaller apertures are more likely to miss 
the fast-moving unbound stars, thereby decreasing the observed 
Mjyn/Mium ratios with respect to larger apertures. 

The outline of the paper is as follows. In Section[2] we quan- 
tify the interplay between the cluster mass function shape and the 
CFRg mass-radius relation in the case of a weak tidal field. Section 
[3] models the time-evolution of the cluster dynamical-to-luminous 
mass ratio, Mj v „/M/„ m . We demonstrate that the time-scale for a 
cluster to recover virial equilibrium after gas expulsion depends on 
the size of the aperture with which it is observed. Our discussion 
helps understand why a few massive clusters are observed in virial 
equilibrium despite a young age of several Myr only. We present 
our conclusions in Section|4] 



1 Note that molecu lar clumps show density gradients ( Miiller et alj|2002l ; 
iBeuther et al]|2002l) . It is thus important to keep in mind that this criterion 
refers to the mean volume density of CFRgs. It does not imply that the 
volume density within a CFRg is uniform. 



2 TRANSIENT CLUSTER MASS FUNCTIONS 

The bound mass of a star cluster after gas expulsion obeys: 

m d =F b x SFE x m C FRg ■ (2) 
© 201? RAS, MNRAS 000,[TL?? 
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Figure 1. Mass-radius diagram of cluster-forming regions (CFRg) with 
lines of constant mean number density («H,.CFRg = 2.3 x 10 4 cm~ 3 , pink 
dotted line), constant mean surface density (XcFRg — 0.5g.crn~ 2 , red solid 
line), constant radius (rcpRg = 0.5 pc, black solid horizontal line), and con- 
stant crossing-times (blue dashed lines with symbols; see key). Illustrated 
are the mass-dependences of CFRg crossing-times in relation with their 
mass-radius relations, hence cluster-mass-dependent evolutionary rates. 



We obtain fi, by linearly int erpolating the /V-body model grid of 
iBaumgardt & Kroupal J2007h rl This provides Fj, as a function of 
t/% cross , SFE, Igexp / Across an d r hl r t (i- e - as m Eq.[T|l. To recover F/,, 
we have to express time in uni t of x cross . For CFRgs with a volume 
densi ty profile of slope —1.7 jMffller et alj|2002l : Parmentier et al. 
l201ll) . we find that the crossing-time at the limiting radius rcFRg 
obeyj^ 



tcross[Myr] = 33 



jrcFRglpc]) 3 
[m C FRg[M e ]) 



1/2 



(3) 



We consider a power-law CFRg mass function of slope —2, 
dAfcffii; 01 m cFR^ m CFR^^ totalizing 10 8 M Q of dense star-forming 
gas. Our aim is to highlight the impact of the sole evolutionary 
rate on the cluster mass function, that is, the impact of t/x cross 
in Eq.[T] To that purpose, we assume all other parameters (SFE, 
r h /r t ) to be mass-independent. We take SFE = 0.35 
and T geX p/t cross = 0.01 (explosive gas expulsion). Note that a 
longer gas-expulsion time-sca le leaves fi, unaffecte d as long as 
Tgexp/tcross IS 0.3 (see fig. 1 in lParmentier et al J2008I) . 

In essence, the tidal field impact r/,/V f is mass-dependent if the 
CFRg mass-radius re lation is not one of constant m ean volume den- 
sity (eqs 8 and 10 in lParmentier & Kroupa ||201 lh . When rcFRg x 

m CFRg ( constant mean surface density), r^/r, <x m^p Rg and higher- 



2 The lib rary of TV-body mode l s on which we build is described in 
detail in IBaumgardt & Kroupa I 120071) . It can be downloaded from: 
http://www.astro.uni-bonn.de/download/data/ The library consists of three 
sets of files: the 'esc'-extension files describe the time evolution of the star 
bound fraction Fj,, the 'rad_bound'-extension files quantify the time evolu- 
tion of the three-dimensional Lagrangian radii of the bound stars, and the 
'rad_bound_2d' -extension files describe the time evolution of the projected 
Lagrangian radii of all cluster stars, i.e. bound and unbound unlike. 

3 Formally, the crossing-time relevant to th e iV-body time-unit is the 
crossing-time at the virial radius of the CFRg (Heggie & Hut 1120031) . It is 
shorter than the crossing-time at the CFRg limiting radius given by Eq.|3] 
by about one third only and we thus neglect that difference. All through 
Sectionf2] the crossing-time of a CFRg is as defined by Eq.fJ] 



Table 1. Adopted mass-radius relations for cluster-forming regions (CFRg). 
Also given are the mean surface density, mean volume density, and radius 
of the T,cFR g , pcFRg and rcpRg models, respectively, as well as the corre- 
sponding CFRg crossing-times, t cross , in dependence of nic FRg . When not 
quoted explicitly, units are pc, Mr., and Myr. 





^CFRg 


PCFRg 


r CFRg 


rcFRg\pc] = 

w[*] = 


0.5g.cm~ 2 
0.04m^ 


1.110~ 19 gxm~ 3 
0.05m^ g 
0.41 


0.5pc 
0.5 

UJm CFRg 



mass clusters experience greater infant weight-loss through their 

—1/3 

tidal boundary. In case of constant radius rcFRg, r h/ r t x m cFRi> an ^ 
higher-mass clusters are more resilient to the external tidal field. To 
'switch off' those mass-dependent effects, we assume that the envi- 
ronment exerts upon star clusters a tidal field weak enough so that 
r h/ r t 0-01 over the entire mass range. When r/Jr, ^ 0.01, clus- 
ters respond to gas expulsion as if in a tidal-field-free environment 
( IBaumgardt & Kr oupa 2007]). In other words, the bound fraction 
F\j does not respond to variations in r/Jr, as long as r^jr, ^ 0.01. 

Figure [TJ illustrates the CFRg mass-radius relations encom- 
passed by our simulations. We refer to them as the ScFSg (con- 
stant mean surface density), PcFRg (constant mean volume density) 
and rcFRt; (constant radius) models. They all lie below the line 
r hl r t = 0.01. As just explained, this cancels any tidally-induced 
mass-dependent effect upon Fj, despite I^cFRg an d r cFR% mod- 
els having mass-varying r/,/r f . Figure [T] also shows lines of iso- 
Tcross ( see Eq. [3}. By virtue of their mass-dependent volume den- 
sities, the rcFRi> an d Ecfr? models give rise to crossing-times 
shorter and longer, respectively, towards higher masses. Model 
properties, along with their corresponding crossing-time in depen- 
dence of mass, are given in Table [T] The adopted mass-radius re- 
lations are representative of star-fo rming molecular clumps (se e 
middle and bottom panels of fig. 1 in lParmentier & Kroupa luOllh . 
We stress, however, that molecular clump sizes are often FWHM 
measurements which, therefore, do not necessarily represent the 
outer boundary of the star clusters formed by these clumps (see 
IParmentierl fcpill) and lParmentier et all ^201 ll) for applications of 
this concept). Other normalizations of the mass-radius relations 
would simply affect T cross hence the age - but not the shapes - 
of the cluster mass functions we obtain below. 

The stellar mass of an embedded cluster obeys: m ec ; = 
SFE x mcFRg- With mass-independent SFE, the embedded-cluster 
mass function is simply the CFRg mass function shifted horizon- 
tally by logio(SFE) ~ —0.5. In addition, with mass-independent 
tgexpl Across, an d fft/Vf < 0.01, the cluster mass function at the 
end of violent relaxation is the embedded-cluster mass function 
shifted horizontally by logio(Fi, ounc i), where Fi, ounc i is the mass- 
independent bound mass fraction at the end of violent relaxation. 
We refer to the bound fraction at the end of violent relaxation as 
Fbound to make it distinct from Fj,, the bound fraction at any time 
during violent relaxation. 

For the model parameters adopted here - SFE = 0.35, 



l gexp I Across 



0-01, r,,/r, < 0.01 - F bound = 0.11. That is, the 
power-law mass functions of CFRgs, embedded clusters and post- 
violent-relaxation clusters all have the same slope. In Fig. [2] they 
are schematically depicted as the (black) dashed lines labelled 
'CFRg', 'eel' and 'VR-end', respectively. Their slope is —1 since 
Fig. [2] represents the number of clusters per constant logarithmic 
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Figure 2. Evolution of the cluster mass function through violent relax- 
ation for three CFRg mass-radius relations: constant mean surface den- 
sity (Y.cFRg)' constant mean volume density (PcFRg) an d constant radius 
U'CFRg), from top to bottom, respectively. Each panel depicts the CFRg 
mass function (solid black line with asterisks) and cluster mass functions 
at three distinct ages (see key). The mass of a cluster is denned as the stellar 
mass enclosed within its tidal radius at the age of relevance. In top and bot- 
tom panels, violent relaxation transiently distorts the cluster mass function 
compared to the CFRg mass function because of cluster-mass-dependent 
evolutionary rates. 

mass bin, i.e. dN/d login. In all this section, the slope of a clus- 
ter mass function refers to the slope of dN/dlogm (e.g. dN c i °= 
mT^Am c i has a slope of — 1 ). 

Over the course of violent relaxation, however, the cluster 
mass function shape is modified if the evolutionary rate is mass- 
dependent. For the £cfr? model, clusters of higher mass evolve 
more slowly than low-mass ones, which renders the cluster mass 
function transiently shallower than the embedded-cluster mass 
function. Conversely, for the r^FRg model, more massive clusters 
evolve faster and we thus expect the opposite behaviour. 

The main results of our simulations are presented in Fig. [2] 



Figure 3. Bound fraction Ft as a function of CFRg mass and cluster age. 
CFRg mass-radius relations and cluster ages are identical to Fig. ~ ] (same 
colour/symbol codings). For the T.cFRg model (top panel), massive clusters 
evolve more slowly than low-mass ones, while the opposite is true for con- 
stant radius (bottom panel). For CFRgs of constant mean volume density 
(middle panel), the cluster evolutionary rate is mass-independent 



with the 1.cfri>, PcFRt> and r CFRg models in top, middle and bottom 
panels, respectively. Each panel presents the 'parent' CFRg mass 
function (solid (black) line with asteriks) and model cluster mass 
functions at three distinct ages (see panel key). To limit the Poisson 
noise of our Monte-Carlo simulations, each displayed mass func- 
tion is the median of a set of 100 simulations. In the top two panels, 
the youngest cluster mass function (red solid line with open circles) 
does not differ from the embedded-cluster mass function. That is, 
clusters are still dynamically too young for the first stars due to be 
unbound to have reached their limiting tidal radius. 

As expected, for the PcFRg model (middle panel), the shape 
of the cluster mass function does not change with time since Z cross 
is mass-independent. For the mean volume density PcFRg assumed 
here, violent relaxation is over by an age of 60 Myr (see also Fig.O, 
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Figure 4. Top panel: Cluster mass function in dependence of age and SFE 
(see key) for CFRgs with a constant mean surface density. In addition to 
lessening the cluster survival rate, a lower SFE is conducive to a greater 
slope change during violent relaxation. Bottom panel: Bound fraction fj in 
dependence of mcFRg for two ages and two SFEs. At an age of 50 Myr, a 
lower SFE stengthens the contrast between the bound fractions at low and 
high masses (compare the blue vertical arrow, SFE = 0.34, to the red one, 
SFE = 0.40), hence the larger slope change at lower SFEs. 



that is, all stars due to be unbound have crossed the cluster in- 
stantaneous tidal radiu.<0. For the CFRg surface density considered 
in the top panel, violent relaxation is completed by clusters of all 
masses by an age of 300 Myr. While violent relaxation proceeds, 
the cluster mass function becomes shallower (slope ~ —0.8; solid 
blue line with plain squares) than the CFRg mass function. The 
opposite behaviour characterizes the r^FRg model (bottom panel), 
with more massive clusters reaching their post-violent-relaxation 
mass at younger ages than low-mass ones. In this case, the transient 
shape of the cluster mass function can be described as a 3-segment 
sequence. Note that the cluster mass functions at the end of vio- 
lent relaxation (green lines with open triangles) are identical for all 
three models. This is because each model is characterized by the 
same CFRg mass function, the same total mass in star-forming gas 
(thus the same CFRg mass function normalization) and the same fi- 
nal bound fraction F\, ounc \ (since all simulations are performed with 
identical SFE, X gexp / 'teross and ri,/r, < 0.01; see Eq.[T}. 

Figure [3] provides another perspective on the same processes. 
It shows the instantaneous bound fraction Fj as a function of mcFRg 
for the same ages as in Fig. [2] In other words, it represents on a lin- 



4 It is worth keeping in mind that the cluster tidal radius itself is shrinking 
with time since it depends on the bound cluster mass, i.e. r t °= m'/ 3 . Hence 
the terminology 'instantaneous tidal radius' adopted here 



ear scale the horizontal shift between each cluster mass function 
and the embedded-cluster mass function as a function of nicFRs;- 
Panel sequence and colour/symbol codings are identical to Fig. [2] 
Figure [3] illustrates both the F/, decrease with time (i.e. as infant- 
weight loss proceeds), and the mass-dependent behaviour of Ej, for 
the 1-cfri> and r CFR<> models. In each panel, the (black) horizon- 
tal dashed line marks the final bound fraction, Ft, olmc i = 0.11 for 
our model parameters. It is reached first by the high-mass clus- 
ters (mass > 10 5 Mq) in the rcFRg model, while the opposite is 
true for the EcFSg model. The mass-dependent evolutionary rate 
is also highlighted at the bottom of each panel, where a time-span 
of 10 Myr since gas expulsion is given in crossing-time units for 
CFRgs of masses 1O 3 M , 10 5 M Q and 1O 7 M . 

The top panel of Fig. [4] illustrates the cluster mass function 
evolution for the Scfr? model with an SFE lower (SFE = 0.34) 
or higher (SFE = 0.40) than in top panel of Fig. [2] At an age of 
50 Myr, SEE = 0.34 (blue line with open triangles) leads to a 
change of the cluster mass function slope greater than SFE = 0.40 
(red line with open squares). The slope is —0.7 for SFE = 0.34, 
while it stays close to — l over most of the cluster mass range 
when SFE = 0.40. This behaviour is explained in the bottom 
panel which shows Fj, as a function of mcFRg- A lower SEE leads 
to a smaller final bound fraction F\, owl ^, thereby enhancing the 
contrast (see the vertical arrows) between high-mass clusters (mass 
> 10 6 Mq), which have just started to give off stars beyond their 
tidal radius) and low-mass ones (mass < 10 4 M,;,, which have 
completed their violent relaxation by an age of 50 Myr). In turn, 
a lower SFE strengthens the change of the cluster mass function 
slope. 

The models leading to the cluster mass functions presented in 
Fig. [2] build on a single variable, namely, the CFRg mass. Other 
model parameters are either fixed (e.g. the cluster age, the CFRg 
star formation efficiency) or depend on the CFRg mass (e.g. the 
CFRg radius for the PcFRg and l^cFRg models). Yet, real star clus- 
ter systems probably form out of CFRgs with a range of SFEs and 
a mass-radius relation noisier than a pure power-law. In addition, 
observed cluster mass functions are integrated over age ranges of 
several Myr to many tens of Myr, depending on the cluster mean 
age. We now investigate how the cluster mass function patterns 
identified in Fig. [2] respond when the simulations include, simul- 
taneously: (i) a cluster age range, (ii) a Gaussian distribution of 
the SFE and f Hi) a power-law CFRg mass-radius relation superim- 
posed with Gaussian noise. We assume a Gaussian SFE distribution 
with a mean SFE = 0.35 and a standard deviation Osf E = 0.05. For 
each CFRg, a Gaussian noise is added to the logarithm of its radius 
predicted by the mass-radius relations used in Fig. [2] We adopt a 
standard deviation of 0.17 in logio(rcFRg)- As a result, > 99% of 
the CFRgs have SFEs between 0.20 and 0.50 (i.e. the 6rj-width of 
the Gaussian distribution), and span a factor of 10 in radius around 
the mass-radius relations of Fig.Q]and TableQ] 

This second set of cluster mass functions is presented as the 
(red) lines with open circles in Fig. [5] The age ranges over which 
they are integrated are: 30-70 Myr, 10-60 Myr and 3-10 Myr for the 
^CFRt;, Pcfs? and rcpRi, models, respectively. For the sake of com- 
parison with the first set of simulations, Fig. [5] also shows the cor- 
responding single-age mass functions (thin solid blue lines) and the 
thick (black) dashed lines 'CFRg', 'eel' and 'VR-end' of Fig.|2] 

Compared to the first set of simulations, the mass function 
features of the rcFRg model are weakened, with the three-segment 
mass function turned into a double-index power-law (compare the 
blue solid line and the red line with open circles in the bottom panel 
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Figure 5. Cluster mass functions (red line with open circles) integrated 
over a given age range (see key) typical of violent relaxation for three 
mean CFRg mass-radius relations: constant mean surface density (< 
T.cFR g >), constant mean volume density (< pcFRg >) and constant radius 
(< rcFRg >), from top to bottom, respectively. Compared to Fig. [5] the sim- 
ulations now add a Gaussian noise to the mass-radius relations of Table [T] 
and include a Gaussian distribution of the SFE [see text for details]. For the 
sake of comparison with the first set of simulations presented in Fig. [2] each 
panel shows the corresponding single-age mass function (solid blue lines) 
and the thick dashed black lines 'CFRg', 'eel' and 'VR-end'. As previously, 
the mass of a cluster is defined as the stellar mass enclosed within its tidal 
radius at the age of relevance. 



of Fig. [5}. For masses lower than 10 3 M Q , the slope of the new 
cluster mass function mirrors that of the CFRgs while, at higher 
masses, it steepens up to — 1 .4. Such a change should still be de- 
tectable in observed cluster data sets. As for the ScFRg model, the 
cluster mass function slopes from the first and second simulation 
sets do not differ significantly (see the solid blue line and the red 
line with open circles in top panel of Fig. [5}. Along the course of 
violent relaxation, the slope of the cluster mass function changes 



by ~ 0.2 (from —1 to —0.8, then back to —1), which is probably 
just at the detectability limit. 

We stress again that none of our simulations include a signifi- 
cant external tidal field. Yet, the tidal field impact may be important 
for the high-mass CFRgs of the LcFRg model. We come back to this 
point later in this section. 

How much slope variation do observed cluster mass func- 
tions show? The evolution with time of the cluster mass func- 
tions in the Solar Neighbourhood and in the L arge Mag- 
ellanic Cloud was inves tigated by IPiskunov et al.l d2008l) and 
Ide Griis & Anders I J20061) . respe ctively. In the Large Magellanic 
Cloud. Ide Griis & Andersl j2006h find that the cluster mass func- 
tion steepens from —0.80 ±0.10 for clusters younger than 30Myr 
to —0.98 ±0.08 for clusters younger than 500 Myr (their fig. 8 
and table 3). A similar t r end - albeit at younger ages - was iden- 
tified by IPiskunov et al.l d2008l) . They report a steepening of the 
Galactic open cluster mass function from —0.66 ±0.14 for ages 
younger than 8 Myr to —1.13 ±0.08 for clusters younger than 
300 Myr (their fig. 6). Although these two examples seem to point 
towards CFRgs of constant mean surface density, they both re- 
qui re strong cautio n ary no tes. The cluster mass estimates inferred 
by IPiskunov et all (2008) are tidal masses, that is, they repre- 
sent the total mass in stars enclosed within the cluster tidal ra- 
dius, r t . As the cluster mass scales with r 3 , any uncertainty on r, 
strongly affects the cluster mass. This may explain why some of 
th eir cluster masses a re higher (3 • 10 s Mq) than the mass estimate 
of lClarketai] ( l2005l) for Westerlund-1, the most massive Galac- 
tic cluster (1O 5 M0). As for the Large Magellanic Cloud, at ages 
younger than 200 Myr, iBaumgardt et alj Jin prep. I) find a cluster 
mass function slope of — 1.32 ±0.1 1, which is significantly steeper 
than that inferred by Ide Griis & Anders I J2006I) (-0.94 ±0.10; see 



IBaumgardt et alJlTn prepl for a full di scussion). Finally, we note 
that lFall. Chandar & Whitmore 1 00091) also find a slight evolution 
of the cluster mass function slope from — 1 . 14 ± 0.03 over the age 
range l-10Myr, to -1.03 ±0.07 over the age range 10-100 Myr. 
This effect is hardly significant compared to the error bars, how- 
ever. 

In addition, we remind the reader that the results presented in 
Figs. [2j5] assume a weak tidal field impact (i.e. r;,/r, < 0.01). The 
impact of a tidal field upon a po pulation of CFRgs is tightly re lated 
to their mass-radius relation (see lParmentier & Kroupal201 ll . for a 
thorough discussion). Let us assume that all CFRgs are exposed to 
the same tidal field (e.g. they are all located in the same region of 
a given galaxy). If the mean surface density of CFRgs is constant, 
more massive objects have a lower volume density and are therefore 
more vulnerable to tidally-induced mass-losses after gas expulsion. 
If the tidal field is strong enough, a deficit in high-mass clusters 
then chara cterizes the end of viol e nt rel axation (see b ottom panel 
of fig. 3 in lParmentier & Kroupa I J201ll) and fig. 4 in lParmentierl 
d2010h ). and the cluster mass function becomes steeper than the 
embedded-cluster mass function. In contrast, CFRgs with a given 
mean volume density all present the same sensitivity to the external 
tidal field, thereby preventing a distorsion of the cluster mass func- 
tion shape through mass-dependent tidally-induced mass-losses. 

The nature of the CFRg mass-radius relation remains heav- 
ily debated. Based on an analysis of the binary population in 
young star clusters, lMarks& Kroupa I J2012h favour a weak mass- 
radius relation (i.e. close to a constant radius). In contrast, 
iFall. Krumholz & Matzner1 ( l2010l) infer a constant-surface-density 
relation based on the SFE - hence the amount of stellar feedback 
- required to cleanse a CFRg from its residual star-forming gas. 
Note that their analysis does not explore the impact of a mass- 
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dependent gas-expulsion time-scale or of a tidal field (see sec- 
tion 4.1 in IParmentier & Kroupa ||201 ll . for a detailed discussion 
of these points). 

We strongly encourage observers to report - without any pre- 
conceived idea - how the cluster mass function evolves with time, 
and to vary the size of the cluster age ranges over which cluster 
mass functions are integrated. This also requests a careful assess- 
ment of the errors affecting the mass function slope and, therefore, 
a good control of the errors affecting individual cluster mass esti- 
mates. 



3 RETURNING TO VIRIAL EQUILIBRIUM: A 
RADIALLY DEPENDENT TIME-SCALE 

3.1 From the cluster central regions to its outskirts: 
understanding the crucial effect of the aperture 

In the previous section, we estimated the time-scale over which 
infant weight-loss is completed in the case of a weak tidal field. 
That is, we assessed how long it takes for all due-to-be unbound 
stars to cross the cluster limiting tidal radius. These time-scales are 
of order several tens of million years. They depend on the CFRg 
crossing-time, and also on the tidal field strength: a stronger tidal 
field is conducive to a smaller tidal radius hence a shorter time-span 
for the unbound stars to cross it. The simulations of the previous 
section being performed under the assumption of a weak tidal field, 
the time-scales derived in Figs l2l5l are upper limits. 

The time-scale needed for the cluster central regions to re- 
turn to virial equilibrium is even shorter. This radial-dependence 
may contribute to explaining why several star clusters are reported 
to be in virial equilibrium despite be ing no older than a few Myr 
(e.g. Westerlund- 1 , hereafter Wd 1; IMengel & Tacconi-Garmanl 
|2009|) . Since star clusters are expected to be out of virial equi- 
librium after gas expulsion, it may be tempting to conclude from 
this type of observations that gas expulsion plays only a minor role 
in the early evolution of star clusters (through e.g. a high SFE, 



e = 0.33 X gexp = X cross r h /r t = 0.03 [BK07 grid] 



adiabatic gas expulsion (t ? , 



» T c , 



s ), or a subvirial state: see 



xpul 

Goodwin I (|20 09|)). In this section, we bui ld on the observations of 



Mengel & Tacconi-Garman 1 d2007l.l2009h to demonstrate that this 



is n ot necessarily the case. 

IMengel & Tacconi-Garman ] l l2009l) estimate the dynamical 
mass of Wd 1, M f / V „, based on the radial velocity dispersion of 
10 massive stars. They find a value comparable to the cluster lu- 
minous mass estimate, M\ um . That is, Wd 1 is observed in virial 
equilibrium, a r esult which may see m unexpected for a 4-Myr 
old star cluster dGennaro et alj [2oTlh if it is born out of a gas 
expulsion process. The discrepancy is lifted, however, when the 
high density of the Wd 1 CFRg - hence its short crossing-time 
- and the limited size of the observed field-of-view are accounted 
for. In what follows, the observed half-light radius, ry,/ ~ 0.86 pc, 
and luminous mass, M lum ~ 10 5 M a , of Wd 1 are taken from 
IMengel & Tacconi-Garman I d2009h. As for the age, we use the 4- 
Myr est imate of Gennaro et alj 1 201 ll) rathe r than the 6-Myr esti- 
mate of Mengel & Tacconi-Garman N2009I) . This equates with a 
time-span of 3 Myr after gas expulsion if gas dispersal took place 
when the cluster was ~ 1 Myr old. Here, it is crucial to keep in mind 
that the return to virial equilibrium of the cluster does not depend on 
the physical time-span since gas expulsion (i.e. the time expressed 
in units of Myr). Rather, it depends on the time-span expressed in 
units of the initial crossing-time, namely, the CFRg crossing-time. 
A cluster can therefore be dynamically old despite a young stellar 
age. 
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Figure 6. Top panel: Evolution of the 3-dimensional Lagrangian radii of 
stars within the instantaneous tidal radius (depicted as the symbol-free 
dotted black line in the bott om panel) for one of the model clusters of 
iBaumgardt & Kroup a 1 2007). SFE, gas expulsion time-scale and tidal field 
impact are quoted on top of the panel. Bottom panel: Evolution of the pro- 
jected (2-dimensional) Lagrangian radii of all stars (bound and unbound). 
In both panels, mass fractions are quoted in the key. Filled rectangles repre- 
sent two evolutionary stages of Westerlund- 1 since gas expulsion (.r-axis), 
and the radius of the field-of-view encompassing the stars observed by 
IMengel & Ta cconi-Garman (2009) (y-axis) (see text for details). 



The 10 stars used to derive the radial velocity dispersion are 
all observed at projected distances less than 2.5 pc from the cluste r 
centre (see both figs. 1 in lMengel & Tacconi-Garman |2007|.|2009|) . 
Assuming a Plummet' profile, the cluster virial radius, r v , follows 
from the projected half-light radius r/,/: r v = 1.3r/, m ~ 1.5 pc, with 
rj lm the 3-dimensional half-mass radius, and r/, m = 1.3r/,/, where 
the factor 1.3 accounts for the projection onto the sky. WithM/, (H1 ~ 
lO 5 ^^), the present-day crossing-time at the virial r adius (eq. 6 
in IBaumgardt & Kroupa 1 12001 llleggie & Hut 1120031) of Wd 1 is 
"cross 

~ 0.25 Myr. As we shall see in Fig. [9] the observed half- 
light radius is aperture-dependent, an effect which complicates fur- 
ther the estimate of the initial crossing-time. 

Figure [6] depicts the time-evolution of the Lagrangian radii of 
a model star cluster with initial parameters: SFE=0.33, r/, /r, = 0.03 
and Xm = 033t cross , w here %m is the gas expulsion time-scale 
as defined by eq. 1 of IBaumgardt & Kroupal d2007l) . Note that 
IBaumgardt & Kroupa I d2007D model gas expulsion as an exponen- 
tial decrease with time of the cluster gas content. Their gas ex- 
pulsion time-scale %m corresponds to the time when a fraction 



5 A si milar value characte rizes NGC 3603 with m c i : 
0.5 pc jRochauetal. l2010t) 



1O 4 M and r lm 
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e~ = 0.37 of the initial gas mass is left. We define T geX p = 3tm so 
that Zgexp corresponds to a residual gas mass fraction of =0.05, 
i.e., the cluster is practically devoid of gas. In Fig. [6] we thus have 
tgexp = Across- The top panel shows the 3-dimensional Lagrangian 
radii of stars within the instantaneous tidal radius. The bottom panel 
represents the projected Lagrangian radii of all stars, bound and 
unbound (i.e. as would be seen by an observer). Units of the x- 
axis, top and bottom y-axes are the initial crossing-time, initial 
half-mass radius and initial half-light radius, respectively. In a first 
step, let us approximate the initial crossing-time and initial half- 
light radius with their observed estimates. The cluster evolutionary 
stage after gas expulsion (3 Myr = V2z cross if T cross = 0.25 Myr), 
and the radius of the field- o f-view containing the stars used by 
iMengel & Tacconi-Garman ] J2009h to derive the radial velocity 
dispersion (2.5 pc/0.86 pc=2.9r/,/, or 2.3r/, m ), are highlighted as 
the left (red) filled rectangles in both panels of Fig. [6] The tar- 
geted region of the cluster is still out of equilibrium Q However, a 
cluster expands and loses mass following gas expulsion. Its initial 
crossing-time is thus shorter than presently observed and the clus- 
ter is now at an accordingly older dynamical stage. Let us make the 
conservative assumption that the cluster has expanded by a factor 2 
only, and let us neglect mass-loss. This implies y- and x-axis units 
smaller by factors 2 and 2\/2, respectively. The rectangle depictin 
the parameter space covered bv lMengel & Tacconi-Garman I J200' 
is thus shifted towards older evolutionary stages by a factor 2.8, 
while its height is increased by a factor 2 (right blue rectangles). 
It now appears that for the inner regions of Wd 1 c overed by the 
field-of-view o f IMengel & Tacconi-Garmanl J2009h . the response 
to gas expulsion may just be over. We analyse in greater detail the 
case of Wd 1 in Sects |3.2l and |3.3l where we demonstrate that the 
hypothesis of virial equilibrium despite a significant gas expulsion 
episode some 3 Myr ago is indeed viable. 

Coeval with the central regions which have returned to equi- 
librium, outer layers keep expanding, thereby bringing more stars 
accross the tidal radius (thick dotted line in bottom panel). The cen- 
tral regions of a cluster may thus be observed in virial equilibrium 
while, on a larger spatial scale, its external layers keep expand- 
ing. This explains straightforwardly why some Galactic massive 
star clusters as young as a few Myr are observed in virial equilib- 
rium (e.g. Wd 1), while clusters with ages 10-60 Myr still show 
a perturbed surface brightness profile at rad ii > 10 pc (e.g M82- 
FNG C1569-A and NGC1705-1; see fig. 1 in lBastian & Goodwin I 
120061) . Owing to their distance, clusters in starburst galaxies are of- 
ten observed with large linear apertures, thereby highlighting the 
shells of unbound stars expanding on spatial scales > 10 pc. In 
contrast, observations of star clusters in the Galactic disc are of- 
ten limited to their few central parsecs (see Sects l3.2 land l3.4l l. We 
now quantify how the size of the aperture affects the observed ra- 
tio of the dynamical and luminous masses of a cluster, M ( / v „/M/„ m . 
In other words, we derive M^ VJ1 /M/„„, as would be done by an ob- 
server. 



3.2 Quantifying the impact of the aperture 

When the aperture is large, the cylinder it defines (i.e. the cylinder 
with the line-of-sight as its main axis) intercepts a greater fraction 



6 Note that the maximum reached by the 10% and 20% projected La- 
grangian radii at times of 15-20 T„- „ (bottom panel of Fig. [6} does not 
imply the absence of radial motions, but the simultaneous contraction of 
the inner regions and expansion of the outer regions 



of the outer shells containing the stars expanding in the vicinity of 
the tidal radius and beyond. As a result, the measured radial ve- 
locity dispersion rises. This simple exercise demonstrates the im- 
portance of the size of the field-of-view: the smaller the aperture, 
the smaller the observed departure from virial equilibrium, and the 
shorter the time-span required for the cluster targeted region to re- 
turn to equilibrium. 

Figure [7j shows the post-gas-expulsion evolution of 
M-dynjMi um . Each panel corresponds to a different set of model pa- 
rameters, i.e. star formation efficiency (e), gas expulsion time-scale 
(tgexp/ Across), tidal field impact (r^/r,) and their corresponding 
final bound fraction {F oom J). The cases covered by Fig. [7] range 
from the complete disruption of the cluster (Ft, olmc i = in top panel 
[a]) to a cluster barely affected by gas expulsion (Fi, olmc i < 1 in 
bottom panel [d]). Note the degeneracy between panels [b] and 
[c]: the lower SFE of panel [b] (SFE=0.33 compared to 0.40 in 
panel [c]) is compensated by its longer gas-expulsion time-scale 
(tgexp = ?cmss compared to T gexp = in panel [c]). Each panel 
considers in turn four aperture radii, r ap = 3.5, 5, 10 and 20 pc, 
which correspond to 4.5, 6.4, 12.8 and 25.6 cluster initial half-mass 
radii ixhmimi — 0.8 pc; see key in top panel). Measuring in our 
simulations the observed half-light radius for apertures smaller 
than 3.5 pc proved inaccurate and results for smaller apertures are 
therefore not presented here. With the aperture and time expressed 
in units of r/, m ,■„,-, and T cross , respectively, the results of Fig. UJcan 
be rescaled to clusters with different initial densities, sizes, and 
apertures. 

The luminous mass M/ um is the total (true) mass contained 
within the aperture. The cluster dynamical mass M^. n is derived 
from the radial velocity dispersion, cjf lK , measured for all the stars 
contained within the aperture. Assuming a Plummer density profile, 
Md yn , a? and r ap are related through: 



_ 32 <Jf os ■ r hp 
yn ~ It G f{t) 



10 



°L ■ r hp t 1 ■ Vl+f 2 



(1+* 2 ) 



2U.5 . 



, (4) 



where r; lp is the projected half-mass radius observed within 
the aperture (equivalent to the half-light radius if light follows 
mass, and to the Plummer radius for a Plummer model), and t = 
r ap/ r hp- For large apertures, i.e. r ao » r;, g , Eq. [4]is reduced to 
the 'canonical' form (see eq. 3 in lFleck et al l2006h : 



M dyn ~ 10- 



7 los ' r hp 



(5) 



As anticipated above, when the aperture is larger, M dyn /Mi um 
is higher and it takes longer for the observed cluster region to re- 
cover Md yn /M{ um ~ 1 (compare the black curves with plus-signs to 
e.g. the pink curves with open diamonds in panels [b-d] of Fig.[7j>. 

The vertical extent of the (blue) shaded area in Fig. |7J illus- 
trates the uncertainty range affectin g the Wd 1 dynamical mass 
estima te, h ence the ratio Mj vn /Mi um . IMengel & Tacconi-Garmanl 
d2007l) and Mengel & Tacconi-Garmanl j2009h estimate Mjy n to 
be 0.63+g^ x 10 5 M© and 1.5+^ x 10 5 M Q , respectively. Be- 
cause the accuracy of the first estimate is hindered by the sm all 
number (4) of stars studied. IMengel & Tacconi-Garmanl 12009J) in- 
creased their sample up to 10 stars. They stress, however, that this 
larger sample includes 5 (presumed) yellow hypergiants whose ra- 
dial pulsations may contribute an increase of the measured veloc- 
ity dispersion. Either estimate thus suffers from its own drawback. 
In what follows, we adopt the average of both values, namely, 
Ll^Oj x 10 5 M?- As for the cluster luminous mass, M\ um ~ 
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Figure 7. Evolution with time of the cluster dynamical-to-luminous mass 
ratio, Mjyu/Miu,,,, for the star s contained in an aperture of ra dius r ap (based 
on the TV-body simulations o flBaumgardt & Kroupal i2007l) ). The dynami- 
cal mass, Miyn, stems from the radial velocity dispersion (Eq.|4), not from 
star proper motion measurements. The size of the aperture is expressed in 
units of the cluster initial half-mass radius, rk m j n « (see key in top panel). 
Gas-expulsion model parameters are quoted at the bottom of each panel, 
along with their final bound fraction Fi, aum i, Note that the gas-expulsion 
time-scale Xgexp is defined as z eexB = 3?m, where •% is the e-folding time 
of gas expulsion of Baumgardt & Kroupa (2007, their eq. 1). The vertical 
line segments correspond to evolutionary stages of Wd 1 based on various 
initial crossing-time estimates and a time-span of 3 Myr since gas expulsion 
(see text for deta ils). The shaded area illustrates the range of uncertainties 
affecting Mjy„ in lMengel & Tacconi-Garman l j2007Ll2009l) 's studies. Note 
that the y-range in panel [a] is much larger than in panels [b-d] . 

10 5 M Q dClark et al.ir2005l ; lMengel & Tacconi-Garman 1120071) , we 
ignore the uncertainties affecting it. The M^/M/,,,,, ratio is thus 
1T+Q5 and is depicted as the (blue) shaded area in Fig. [7] The 
vertical (black) line segments to the left of panels [b-d] indicate 
the evolutionary stage of Wd 1 after gas expulsion if the initial 
and observed crossing-times are equal, i.e. as for the red (left- 



Figure 8. Evolution of the ratio between the visual luminosity, L,., of a star 
cluster and its dynamical mass, Mj v „, with both quantities measured over 
a cylindrical aperture of radius r ap . These diagrams are the counterparts of 
those in Fig. [7] (identical colour/symbol codings). That is, the time evolu- 
tion oiMjy„/Mi um has been combined with the age-d ependent mass-to -light 
ratio of a single age stellar population according tojMaraston (2()o3). We 
assume that gas is expelled at an age of 1 Myr and that the initial crossing- 
time is 0.15 Myr. Note that time in this figure is measured in years, while 
the time-unit in Fig.UJis the cluster initial crossing-time. 



ward) rectangles in Fig.[6] t = 3Myr = 12x crojs . The line segment 
thick part indicates the difference between the observations (i.e. the 
horizontal dotted line at M^ v „/M; um = 1.1) and the models with 
r ap /rhm,mit=4.5 and 6.4 (see Sect.|33J. 

The amount of dynamical evolution experienced by a cluster 
after gas expulsion depends on its initial crossing-time, namely, the 
crossing-time of the CFRg out of which it formed and defined by 
Eq. [3] The observed crossing-time of a gas-free cluster is longer 
than that of its CFRg because of (i) the gas loss, (ii) the infant 
weight-loss, and (Hi) the gas-expulsion-driven cluster expansion. 
In other words, using the observed crossing-time leads us to under- 
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Figure 9. Evolution of the (projected) half-light radius observed within 
different apertures for the quoted SFE £, gas-expulsion time-scale 1 geX p 
and tidal field impact ri,/r,. The aperture sizes in units of the initial 
half-mass radius, r a p/r}mmit< are given in the key. In the left part of the 
panel, the horizontal line segments indicate the three smallest apertures 
in units of the initial half-light radius, r a p/fhl,tnif F° r the N-body model 
of iBaumgardt & Kroupa 1 120071) under scrutiny here, ruinit — 0.6 pc and 
r ap =3.5, 5, 10 and 20 pc which gives r ap /r^ injt =5.9, 8.5, 17 and 34. The 
last and largest aperture is beyond the plot top border. 



estimate the cluster dynamical age, an effect introduced in Fig. [6] 
In addition, a consistent comparison between the models and the 
observations requires the size of the aperture in units of the initial 
half-mass radius, r ap /r nm j„i t , which is not accessible to observa- 
tions. We now derive estimates of the initial crossing-time and ini- 
tial half-mass radius. 



3.3 An attempt to recover the initial crossing-time 

The 



ra dius of the field o f 
iMengel & Tacconi-Garman I 



iMengel & Tacconi-Garman 



view of Wd 1 

d2007h and 
~ 2.5 pc. Since the 
initial half-mass radius is unknown, we resort to rough esti- 
mates of the ratio r ap /ri lm j„j t in the first place. We assume 

hm,init ~ 2.5. In 

4.5 and 6.4. 



0-3 pc < Own < 1 pc which leads to 8 > r ap /r/. 
Fig-El we therefore consider models with r ap jr nm i n i t 

To answer the question as to whether Wd 1 can be observed in 
virial equilibrium despite a major gas-expulsion phase some 3 Myr 
ago, we focus on panels [b] and [c] of Fig. [JJ where the final in- 
fant weight-loss is significant, i.e. two-thirds of the cluster initial 
stellar mass are lost. (Panel [a] presents the case of a fully dis- 
rupted cluster, irrelevant here, and we come back to panel [d] later 
in the discussion). Building on the observed crossing-time, we find 
a discrepancy between the models (green curve with asterisks and 
pink curve with diamonds) and the observations (horizontal black 
dotted line) larger than 3a, as quantified by the thick part of the 
most leftward vertical segment. As we saw in Fig. [6] however, us- 
ing the observed crossing-time underestimates the cluster evolu- 
tionary stage. So let us now correct the observed crossing-time for 
the gas-loss (by multiplying it by SFE 1 / 2 ; see Eq. [3}- This leads 
to the cluster evolutionary stages indicated by the second vertical 
segment from left in panels [b-c] at t ~ 20z cross . This sole correc- 
tion brings the gas-expulsion models within 1 rr of the observations 



when r ap < lr ht 



In addition to gas loss, the cluster is more 



expulsion time-scale and tidal field impact (same parameters as in 
panel [c] of Fig. [JJ and model outputs similar to panel [b]). The 
observed half-light radius is - in essence - aperture-limited. The 
apertures used to observe the model cluster are shown as horizon- 
tal lines in the left part of the plot with identical colour/symbol 
codings. Note that they have been rescaled in units of the ini- 
tial half-light radius and that the largest aperture is out of range. 
For the apertures and dynamical evolutionary stages of relevance 
here (8 > r ap / ! rhm,init 2.5, t > 20T crosj after correcting for the 
gas loss), the cluster has expanded by a factor of ~ 2. The ini- 
tial half-light radius m ay thus have been >-;,; ,-,„> ~ 0.43 pc (half 
the value measured by IMengel & Tacconi-Garman~ll2007l) which, 
for a Plummer profile, corresponds to an initial half-mass radius 
r hm,init — 0.56 pc (since ri m ini , = 1. 3 r/,/ ,„,,). The factor of two in 
expansion suggests that the initial crossing-time may have been as 
short as r cross ~ 0.25 Myr x SFE 1 / 2 x (1/2) 3 / 2 ~ 0.06 Myr. The 
corresponding dynamical ages after gas expulsion are shown as 
the vertical line segments at t/t cross ~ 55 — 60 in panels b-c of 
Fig. UJ The observed M f / vn /M/„„, ratio and inferred evolutionary 
stages t/tcross of Wd 1 agree neatly with a gas expulsion sce- 
nario in which the stellar mass of the cluster after revirialization 
amounts to one-tenth only of the CFRg mass (since e ~ 0.33 and 
Fbound — 0.33). The ~ 70%- infant weight-loss is visible as the de- 
crease with time of the cluster tidal radius in the bottom panel of 
Fig. [6] r, oc (ffl c( )V 3 oc (Ff, x m ec i) 1 '^, or an eventual decrease of 
0.30 1 ' 3 . It is therefore worth emphasizing that star clus- 



,4/3 
! bound 



compact prior to gas expulsion. Figure [9] illustrates the evolution 
of the observed projected half-light radius for the quoted SFE, gas- 



ters with stellar ages of several Myr, which formed out of dense 
molecular clumps with short crossing-times, are already dynam- 
ically evolved, that is, the time-span since gas expulsion can be 
many initial crossing-times. 

The value inferred for the initial crossing-time at the virial 
radius of Wd 1, r cross ~ 0.06 Myr, may seem short. Still, it is 
not too dissimilar from th e crossing-time o f the m olecular clump 
GO.253+0.016 for which lLongmore etal] J2012h find x cross ~ 
0.17 Myr (their table 2). We note that the mass of GO.253+0.016 
and Wd 1 are similar (i.e. ~ 10 5 Mq). If gas-expulsion played a 
major role in the early evolution of Wd 1, then Wd 1 formed out of 
a CFRg even more massive than GO.253+0.016. With e ~ 0.33 and 
Fbound — 0.33, the parent CFRg of Wd 1 may have been as massive 
as 10 6 M o . A crossing-time of T cross ~ 0.06 Myr is then achieved 
if r^FRg — 1 -5 pc (see Eq.[3]and Fig.Ql. This is r eminiscent of th e 
radius put forward bv lKroupa & Boilv U2002I) and lKroupaU2005l) . 
who assume that the CFRg radius is about 1 pc regardless of the 
CFRg mass. We now see a connection between the time-evolution 
of the cluster mass function shape (Sect. |2j and how quickly star 
clusters return to virial equilibrium, the link between both aspects 
being the CFRg mass-radius relation. 

What if the mass-radius relation of CFRgs is not one of 
constant radius? If CFRgs have a const ant mean volume density 
jWu et al J [20051 : IParmentier et al.l 1201 ll) . then, at high mass, the 
field-of-view may not be large enough to cover the CFRg spatial ex- 
tent in full. For instance, if the mean volume density is that adopted 
in Fig.Q](i.e. n H2 ~ 2.3 • 10 4 cm" 3 ), a CFRg 1O 6 M in mass now 
has a radius rcFRg — 5 pc and a crossing-time % C ross — 0-4 Myr. The 
dynamical evolution is thus markedly slower than the case detailed 
above for rcFRg — 1.5 pc as a time-span of 3 Myr after gas expul- 
sion now corresponds to ~ 8 V cross only. Yet, that does not neces- 
sarily imply that the cluster cannot be observed as having returned 
to equilibrium. The aperture used to observe Wd 1 (~ 2.5 pc) is 
now twice as small as the size of our putative CFRg. This limits 
the observations to the inner half of the parent CFRg, with two im- 
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mediate consequences. Firstly, since molecular clumps have den- 
sity gradients, the initial crossing-time of the observed field-of- 
view is shorter than the 'global' value (Across — 0.4 Myr) quoted 
above. Secondly, assuming a slope of —1. 7 to —2 for the po wer- 
law density gradients of molecular clumps dMuller et alj2002h . wc 
get r ap ~ ronton (that is, the aperture corresponds to about half of 
the three-dimensional CFRg mass), a case not investigated in Fig.[7j 
New models where the aperture covers a fraction only of the CFRg 
are therefore needed, a study which we defer to a future paper. 

Another possibility is of course that gas expulsion is a minor 
perturbation to the cluster and the final bound fraction Ft, ounc i is 
close to unity. This is the case for instance because gas was ex- 
pelled on an adiabatic time-scale as in panel [d] of Fig. [7] The two 
vertical line segments depict the evolutionary stages of Wd 1 based 
on the observed and gas-loss-corrected crossing-times. In case of 
adiabatic gas expulsion, the embedded cluster experience neither 
significant infant weight-loss nor significant spatial expansion and 
further corrections of the observed crossing-time are therefore su- 
perfluous. The limited cluster expansion allows the cluster to return 
to equilibrium on a time-scale as short as — &T C i- oss for apertures 
smaller than 10r/ lm ,,„y (and about 3T cm ss at the lcr level). 

To estimate how quickly the intra-cluster gas is removed from 
an embedded cluster is a complex problem. If the ionized mate- 
rial simply expands at the sound speed, vs ~ 10km ■ s , then a 
CFRg of mass mcFRg — 10 6 Mq and radius ''CFRg — 1-5 pc loses 
its gas on a time-scale of 0.15 Myr, equivalent to a few crossing- 
times (sin ce T cross ~ 0.06 Myr; s ee above). Gas expulsion is thus 
adiabatic dKroupa &Boilvll2002l) and the return to equilibrium is 
best-described by panel [d]. The gas-expulsion time-scale can be 
shorter {x gex p < T cross ), ho wever, if radiation pressure d ominates 
the gas-expulsion process dKrumholz & Matzner I [2009h and the 
young cluster dynamical evolution then follows from panels [b-c]. 

In Fig. [7] the vertical extent of the shaded (blue) area ig- 
nores the uncertainties affecting M\ um . Were they included, the 
error on M ( / v „/M/„ m would be accordingly larger. In that respect, 
iGennaro et ail feOl lh derive M/, ,,„ ~ 0.49+° ^ x 10 s M @ for Wd 1 
(a similar value was fou nd by iBrandner et al. 1120081) . Com bined 
to the Mctm estimate of iMengel & Tacconi-Garman I d2009l) . this 
is suggestive of a higher ratio Mj m / 'Mi um ~ 2 (i.e. in Fig. [Vj the 
shaded areas rise by a factor of ~ 2), even more compatible with 
explosive gas expulsion conditions. Here lies one more source of 
complication, however. In models, the stars used to derive the ra- 
dial velocity dispersion a? and the luminous mass estimate Mi um 
are identical. In real observations, that may not be the case. To infer 
M[ wn and M$y n based on different observation sets, hence different 
fields-of-view, can only increase the error in M,; v „/M;„ m . 

Figure [8] combines the time-evolution of the M c j xn /Mi um ra- 
tio of Fig. [7] with the age-dependent mass-to-light rati o of a Sim- 
ple Stellar Population as predicted bv lMaraston I |2Q05b for a solar 
metallicity. We assume that gas expulsion occurs when the cluster 
is 1 Myr old. Because the dynamical and photometric evolutions 
of star clusters have different time-units (initial crossing-time and 
Myr, respectively), Fig. ^demands a reliable estimate of the ini- 
tial crossing-tim e. We take z cros s = 0. 15 Myr, a value similar to 
that inferred by lLongmore et al .1 d2012l) for the molecular clump 
GO. 253+0.016. T cross = 0.15 Myr is also intermediate between our 
estimate of the initial crossing-time of Wd 1 and its observed value 
(0.06 Myr and 0.25 Myr, respectively). Panels, shaded areas and 
evolutionary tracks in Fig.[8]are the counterparts of those in Fig. [7] 
(identical colour/symbol codings). 

Some models of the ti me-evolution of M r j y „ /M/u m and 
Ly/Mjyn have been presented in lGoodwin & Bastian U2006I) . Their 




Figure 10. Mapping the gas-expulsion driven expansion of a cluster. The 
line-of-sight is indicated by the (black) arrow on the left. Radial and proper 
motions are depicted by the horizontal (blue) and vertical (red) arrows. The 
aperture is shown as the shaded area. Because of the narrow centrally- 
located aperture, a proper motion mapping within the aperture misses the 
expansion of the cluster outer layers. In contrast, radial velocities trace the 
whole radial extent of the cluster, from the central regions to the outskirts. 
As a result, the time-evolutions of the dynamical masses inside the aper- 
ture inferred by proper motions on the one hand and radial velocities on the 
other hand are expected to differ. 



simulations build on a Plummer model with a Plummer radius 
rp = 3.5 pc and a star+gas mass mcFRg — 5 x 1O 4 M0. This yields 
a virial radius of r v ~ 1 , 7 ■ ~ 5 . 8 pc and a crossin g-time T C ross = 
2.6 Myr (using eq. 6 in lBaumgardt & Kroupa 1120071) . Compared to 
the initial crossing-time of a starburst cluster, this is about an or- 
der of magnitude too long. We also note that their dynamical mass 
is measured from the velocity dispersion of stars located within 
20 pc of the cluster centre (equivalent to 4.5r/ lm lmf in their model). 
In contrast, our velocity dispersion builds on the stars contained 
within a cyli ndrical aperture. The r efore, given that (i) it is un- 
clear whether iGoodwin & Bastian 1 d2006l) rescaled their model to 
a crossing-time appropriate for starburst clusters, and ( ii) we define 
the M e i vn /Mi lm ratio based on different ensembles of stars, our re- 
sults cannot be consistently compared to theirs. Additionally, the 
tidal field strength is set to zero in their simulations, while Figs l7l8l 
builds on a weak tidal field. 



3.4 The Arches cluster: proper motions and extreme 
aperture 

A recent study of the Arches cluster by Iciarkson et al .1 d2012l) in- 
fers estimates for its dynamical and luminous masses which are 
similar. B ecause the Arches c luster is younger than Wd 1 (age: 2- 

2.5 Myr, iNaiarro et alj|2004h . the time-span since gas expulsion 
may be as short as 10 initial crossing-times (assuming an age of 
1 Myr at gas expulsion and % C mss — 0.1 Myr). At first glance, this 
seems to disprove a gas-expulsion scenario like that shown in pan- 
els b-c of Fig. |7j However, it proves crucial to consider how the 
dyn amical mass is derived and, again, how small the aperture is. 

Iciarkson et al. I d2012h map proper motions of stars within 
the central 10" x 10"of the cluster (or 0.4 pc x 0.4 pc, equivalent 
to a projected radius of 0.2 pc). Combining their proper motion 
measurements with cluster modelling, they estimate the dynami- 
cal mass contained within a cylinder of radius 0.4 pc to be about 
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0.9 • 10 4 Af Q . This is comparable to their luminous mass estimat^ 
within the same aperture, namely, M\ llm =1-2-10 4 Mq. The small 
radius of 0.4 pc is comparable to, or smaller than, the initial half- 
mass radius of the cluster (see our estimate for Wd 1 above). One 
thus needs a model with r ap /r^ m j n i t < 1, a limit which our present 
,/V-body simulations do not achieve. All we can anticipate is that, 
within such a small aperture, the cluster is observed to return to 
virial equilibrium even faster than in any of our simulations. There 
is a strong caveat, however. In the models presented in Fig. [7] 
the dynamical mass is estimated based on radial velocity mea- 
surem ents (see Eq. [4jl. The Arches cluster study of Iciarkson et al.l 
42012b builds on proper motion measurements. Although the devel- 
opment of a model akin to Eq. [4] for proper-motion-based dynam- 
ical masses is beyond the scope of the present paper, we point out 
that the measurements of proper motions in a small aperture neces- 
sarily miss the expanding external layers which characterize a star 
cluster after gas expulsion. This is illustrated in Fig. [10] where the 
horizontal (blue) arrows depict the radial motions and the vertical 
(red) ones depict the proper motions. The aperture is shown as the 
shaded area. While radial velocity measurements probe the full ra- 
dial extent of the expanding cluster, proper motion measurements 
are restricted to the 3-dimensional central regions. One may thus 
expect a lower measured velocity dispersion when based on proper 
motions than when based on radial velocities, hence M% < 
where M^'" n and are the dynamical mass estimates obtained 
from the proper motions and radial velocities, respectively. To sum 
up, because of an aperture comparable in size to the initial half- 
mass radius of the cluster, and because of a cluster dynamical mass 
estimated from the proper motions of the cluster central regions 
only, it cannot be excluded that M^™ ~ M; Hm is observed at the 
present Arches-cluster age even though the cluster experienced a 
major gas expulsion episode in the recent past. 

3.5 What must be kept on top of mind 

The above discussions highlight the complexity of inferring 
whether the observed virial equilibrium of young gas-free star 
clusters is compatible with a gas-expulsion scenario. A key-issue 
is that pivotal model parameters such as (1) the initial crossing- 
time, hence the CFRg mean volume density, and (2) the size of 
the field-of-view with respect to the cluster initial half-mass radius, 
are not well-constrained. It is of course crucial to ensure that mod- 
els and observations are comparable. That is, a model with a long 
cros sing-time might be appropriate for the so-called 'leaky' clus- 
ters jPfalznerll2009h , while it severely underestimates the rate of 
dynamical evolution of 'starburst' clusters such as Wd 1 and the 
Arches. 

The observed M^ n jmi um ratio is subject to its own un- 
certainties, related to (1) the cluster luminous mass estimate, 
Mium ( see the exemple of Wd 1 in the above discussion), (2) 
the measured velocity dispersion, of , and (3) the observed 
half-light radius. Here, an additional source of complication is 
that the half-light radius is stellar-mass dependent if the clus- 
ter is primordial ly mass-segregated ( which may be the case for 
NGC 3603; see (Rochau et adhoicl) . Note that the models pre- 
sented here are not primordially mass-segregated since the N- 

1 Their photometric mass estimate is derived under the assumption of a 
top-heavy initial stellar m ass function, an issue which remains disputed for 
the Arches cluster (see e.g. lEspinoza. Selman & M elnick 2009j), and which 
is beyond the scope of the present discussion. 



body simulations of iBaumgardt & Kro upa (2007j) build on equal- 
mass particles. [Models for the evolution of primordially mass- 
segregated gas-free clusters are presented in iFleck et al 
and the impact of stellar-evolutionary mass- losses is studied in 
IVesperini. McMillan & Portegies Zwartl J2009h l. 

Figure [7j clearly demonstrates that star clusters as young as a 
few Myr observed in virial equilibrium and gas-expulsion driven 
cluster evolution are not mutually exclusive. We encourage ob- 
servers to estimate and quote the crossing-time of the star clusters 
they study, as well as the size of the field of view they scrutinize. 
Both parameters play a pivotal role in assessing reliably the dynam- 
ical state of a star cluster. Finally, we caution that in Figs [7] and [8] 
the cluster dynamical mass is obtained from the radial velocity dis- 
persion. Models of the time-evolution of /Mi um , where M^'" n 
builds on star proper motions, are highly desirable. 



4 FUTURE WORK AND CONCLUSIONS 

We have studied the impact of cluster-mass-dependent evolutionary 
rates upon the evolution of the cluster mass function through vio- 
lent relaxation. To highlight this so far overlooked process, we have 
built model star cluster systems with mass-independent star forma- 
tion efficiency (SFE), gas-expulsion time-scale (jgexp/ Across) and 
tidal field impact (r/,/r f ). Due to the resulting mass-independent 
infant weight-loss at the end of violent relaxation, the final clus- 
ter mass function and the embedded-cluster mass function have 
the same shape. Yet, that does not necessarily imply that the clus- 
ter mass function during violent relaxation retains the same shape 
as that of embedded clusters. If the evolutionary rate is mass- 
dependent, it transiently distorts the cluster mass function. 

Mass-dependent evolutionary rates arise when the CFRg 

mass-radius relation is not one of constant mean volume density 

— 1/2 

(since T cross oc p CFR ). For instance, for CFRgs with a constant 
radius (rcFRg model, bottom panels of Figs I2l3t , more massive 
clusters evolve faster and complete their violent relaxation earlier 
than low-mass clusters. Conversely, in case of CFRgs of constant 
mean surface density, dynamical evolution is slower in the high- 
mass regime (X-CFRg model, top panels of Figs [2~T3b . Therefore, to 
preserve the shape of the cluster mass function during violent relax- 
ation requires CFRgs of constant mean volume density. In that case, 
the bound fraction of stars is mass-independent at any age during 
violent relaxation (PcFRg model, middle panels of Figs I2l3t . We 
caution, however, that the influence of the CFRg mass-radius re- 
lation upon the evolving cluster mass function is weakened when 
the scatter in cluster initial properties is accounted for (e.g. range 
in SFE, scatter in radius around a mean mass-radius relation, etc; 
see Fig. [5}. 

A direct comparison of our models to observations still re- 
quires additional work, however. In our simulations, the mass of 
a cluster is defined as its bound mass, namely, the mass enclosed 
within its instantaneous tidal radius. This may differ from the 
cluster mass inferred by observers. Due to their gas-expulsion- 
driven spatial expan sion, star clusters develop untruncated power- 
law density profiles telsonetalJI 19871) at young ages. That is, un- 
bound stars beyond the tidal radius are not yet spatially dissoci- 
ated from their 'parent' clusters. Therefore, contrary to older clus- 
ters, no truncation in the cluster density profile marks its tidal ra- 
dius. The observed mass of those young clusters is thus inferred 
as the clus ter stellar mass bright er than the host galaxy back- 
ground (see lLtighausen et alj2012l . for an introductory work). Nev- 
ertheless, infant weight-loss beyond the tidal radius and cluster 
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spatial expansion are correla ted (see e.g. iGever & Burkert 11200 ll ; 
iBaumgardt & Kroupa Il2007h . At a given age, stronger spatial ex- 
pansion dims the cluster surface brightness. Therefore, the cluster 
bound mass fraction and the cluster mass fraction brighter than the 
background of stars and gas of the host galaxy are correlated too. 
We thus expect the same cluster mass function patterns for cluster 
masses defined based on a surface brightness limit criterion, al- 
though the amplitude of the mass function distortions may differ 
from those found here. In spite of this caveat, it is clear that our 
work has highlighted, again, the richness of the cluster mass func- 
tion as a tracer of star cluster formation conditions. 

On the scale of individual clusters, violent relaxation can be 
traced by the evolution of the ratio between the dynamical and lu- 
minous masses of a cluster, Mj xn /Mi um . Models of M^ yn /Mi um as 
a function of dynamical time have been computed for various star 
formation efficiencies and gas-expulsion time-scales (Fig. [7}. Our 
simulations range from disrupted clusters to clusters barely affected 
by gas expulsion. We insist that for a given set of parameters (SFE, 
tgexp I Across, ^h/ r t)< me rate of dynamical evolution is determined 
by the CFRg crossing-time. 

Our study has quantified an effect ignored so far: the impact 
of the size of the field-of-view on the observed M,; v „/M/ Mm ra- 
tio. Smaller apertures are conducive to smaller Af</ yn /M/ Hm ratios 
for the following two reasons (1) Cluster inner regions return to 
equilibrium at a time when outer layers still expand through the 
tidal radius (Fig. |6j. (2) A smaller aperture 'sees' a smaller frac- 
tion of the stars expanding beyond the cluster tidal radius (i.e. at 
large distance from the cluster centre, the only expanding stars to 
be observed are those in the vicinity of the line-of-sight), which 
decreases the measured radial velocity dispersion a 2 . Both effects 
lower the mass ratio Md yn /Mi um , that is, a cluster observed after 
gas expulsion through a smaller aperture is perceived as closer to 
dynamical equilibrium and its return to equilibrium is seen as faster 
(Fig. [7). We stress that what actually matters is the size of the aper- 
ture in units of the cluster initial half-mass radius. If the cluster 
initial crossing-time is short and the observed field-of-view is no 
larger than a few pc, it is indeed possible to infer M^ yn ~ M[ wn 
even though the cluster experienced significant gas expulsion a few 
Myr earlier. 

Both the cluster mass function and cluster dynamical mass as- 
pects of our work emphasize the crucial importance of the crossing- 
time to the evolutionary rate of young star clusters. We strongly 
encourage observers to estimate and quote the crossing-time of the 
young star clusters they study, although the observed crossing-time 
may differ from the initial one by a factor of a few. From a dynam- 
ical point of view, the age of a star cluster expressed in physical- 
time units hardly tells anything about its evolutionary stage. The 
relevant quantity is the age expressed in units of the cluster initial 
crossing-time, that is, the CFRg crossing-time (see Fig. |T}. To ob- 
tain a reliable estimate of the latter is thus crucial when it comes to 
combining dynamical and photometric modellings of star clusters, 
as their respective time-units, the initial crossing-time and the year, 
differ (e.g. see Fig.[8]for a conversion of the M,; v „/M/„ m ratio into 
the luminosity-to-dynamical mass ratio). 

The mass-radius relation of CFRgs - both its slope and nor- 
malization - determines the initial crossing-time of star clusters. As 
such, it remains central to our understanding of star cluster dynam- 
ical evolution after residual star-forming gas expulsion. 
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